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ON  THE  STABILITY  OF  THE  LAMINAR  MIXING  REGION 
BETWEEN  TWO  PARALLEL  STREAMS  IN  A  GAS 
By  C.  C.  Lin 


SUMMARY 


The  stability  of  the  mixing  of  two  streams  was  studied  both  for 
the  interest  in  the  problem  and  for  clarifying  certain  points  in  the 
basic  stability  theory.  It  is  shown  that,  when  the  relative  speed  of 
the  two  parallel  streams  exceeds  the  sum  of  their  velocities  of  sound, 
subsonic  oscillations  cannot  occur  and  the  mixing  region  may  be  expected 
to  be  stable  with  respect  to  small  disturbances.  It  is  further  shown 
that,  when  viscosity  and  heat  conductivity  are  neglected,  if  the  flow 
can  execute  a  small  neutral  subsonic  oscillation  of  finite  wave  length, 
it  can  also  execute  self -excited  oscillations  of  longer  wave  lengths 
and  damped  oscillations  of  shorter  wave  lengths. 

Rigorous  developments  of  the  mathematical  theory  of  asymptotic 
solutions  confirm  previous  methods  of  solution  of  the  stability  equa¬ 
tions  in  a  compressible  fluid.  This  theory  also  shows  that,  at  high 
Reynolds  numbers,  the  damped  oscillations  in  a  strictly  parallel  main 
flow  have  a  structure  similar  to  that  of  the  vorticity  field  in  fully 
developed  turbulent  flow. 

Sample  calculations  are  also  included  exhibiting  various  quantita¬ 
tive  properties  of  these  small  oscillations. 


INTRODUCTION 


The  mixing  of  two  parallel  streams  of  gas  occurs  in  a  number  of 
cases.  An  interesting  example  Is  furnished  by  the  slip  stream  in  a 
three-shock  configuration.  It  has  long  been  suggested  that  such  laminar 
mixing  zones  could,  at  sufficiently  high  speeds,  be  stable  with  respect 
to  small  disturbances  although  they  are  known  to  be  very  unstable  at  low 
speeds.  The  purpose  of  the  present  investigation  is  to  find  out  some 
of  these  stability  characteristics.  Apart  from  the  development  of  the 
general  theory,  there  are  included  the  calculations  of  the  neutral  and 
unstable  oscillations,  the  extent  of  the  amplification,  and  other  related 
properties. 
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The  basic  equations  for  the  study  of  small  disturbances  in  the 
laminar  boundary  layer  of  gases  have  been  given  in  reference  1.  However, 
the  development  of  the  theory  there  has  specific  reference  to  the  case 
of  a  layer  near  a  solid  boundary.  As  pointed  out  in  references  2  and  3, 
the  stability  theory  for  a  mixing  zone  in  an  incompressible  fluid  dif¬ 
fers  from  that  for  a  layer  near  a  solid  boundary  in  that  solutions  of 
the  exponential  type  are  unimportant.  This  leads  to  the  conclusion  that 
the  effects  of  viscosity  and  heat  conduction  are  negligible,  except  at 
very  low  Reynolds  numbers,  in  determining  the  characteristics  of  the 
oscillations.  To  confirm  this  point,  a  rigorous  mathematical  theory  of 
asymptotic  solutions  was  developed  for  the  compressible  case  similar  to 
that  indicated  in  reference  4  for  the  incompressible  case.  In  view  of 
the  mathematical  interest  involved,  it  was  decided  that  this  basic  part 
of  the  present  investigation  would  be  published  separately  in  mathematical 
journals  (see  references  5  and  6),  and  only  the  main  results  and  their 
physical  significance  will  be  presented  here. 

As  in  the  case  of  an  incompressible  fluid,  the  "inviscid"  case  is 
expected  to  be  characteristic  of  the  behavior  of  the  disturbances  at 
moderately  large  Reynolds  numbers.  Mest  of  the  studies  are,  therefore, 
made  in  the  inviscid  case.  However,  the  interpretation  of  the  inviscid 
case  must  be  subjected  to  the  same  care  as  in  the  incompressible  case; 
that  is,  in  the  case  of  damped  disturbances,  the  differential  equation 
of  the  inviscid  flow  may  not  be  regarded  as  valid  throughout  the  real 
axis.  There  is  a  finite  viscous  region  even  in  the  limit  of  vanishing 
viscosity.  The  complex  conjugate  of  the  amplified  disturbance  is  cer¬ 
tainly  a  solution  of  the  inviscid  equation,  but  it  is  not  a  limiting 
solution  of  the  complete  viscous  equation.  This  behavior  of  the  inviscid 
solution  reminds  one  of  the  vorticity  structure  of  fully  developed  tur¬ 
bulent  flow  as  found  by  Batchelor  and  Townsend  (reference  7) • 

The  nonexistence  of  subsonic  disturbances  is  visually  associated 
with  the  stability  of  the  parallel  flow.  There  seems  to  be  some  basis 
for  doing  this,  although  the  role  of  supersonic  disturbances  has  never 
been  fully  clarified.  It  is  easy  to  see  that,  for  certain  combinations 
of  the  properties  of  the  two  streams,  it  is  impossible  to  have  a  sub¬ 
sonic  disturbance  relative  to  both.  Under  such  conditions,  one  may 
expect  stability.  These  conditions  for  stability  are  developed  herein 
and  are  expected  to  hold,  irrespective  of  the  viscous  effects. 

Applying  the  theory  of  stability  in  the  inviscid  case,  one  can 
further  narrow  down  the  possible  range  of  instability.  This  will  depend 
upon  the  velocity  and  temperature  distributions  in  the  shear  zone.  In 
the  present  work,  calculations  are  made  for  gases  with  Prandtl  number 
equal  to  unity.  Although  the  condition  of  equal  total  enthalpy  in  the 
two  streams  is  also  used,  it  is  shown  that  this  restriction  can  be 
immediately  removed  by  considering  a  moving  observer.  It  is  found 
that  the  condition  of  stability  thus  found  does  not  differ  very  much 
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,  from  that  found  above  from  general  considerations.  Thus,  it  may  be 

surmised  that  the  exact  distributions  of  temperature  and  velocity  have 
only  a  secondary  influence  on  the  stability  characteristics  in  the 
*  mixing  zone.  Thus,  the  approximations  used  in  the  present  calculations 

of  the  basic  velocity  and  temperature  distributions  cannot  influence 
the  final  results  to  any  appreciable  extent. 

The  neutral  disturbances  are  of  two  kinds:  (l)  A  steady  deviation 
and  (2)  an  oscillation  of  finite  wave  length.  Thus,  there  are  two 
branches  of  the  neutral  curve  at  infinite  Reynolds  numbers.  They  may 
be  expected  to  join  together  at  low  Reynolds  numbers  enclosing  a  region 
of  instability.  Calculations  of  neutral  and  amplified  oscillations  are 
carried  out  in  a  number  of  cases  with  one  stream  at  rest. 

This  investigation,  carried  out  at  the  Nkssachusetts  Institute  of 
Technology,  was  sponsored  by  and  conducted  with  the  financial  assistance 
of  the  Rational  Advisory  Committee  for  Aeronautics.  The  author  is 
indebted  to  Mr.  D.  W.  Dunn  for  his  valuable  suggestions  and  help  in  the 
preparation  of  the  final  version  of  the  report  and  to  Miss  Diana  Ifeson 
and  Mr.  W.  V.  Caldwell  for  their  help  in  making  the  numerical  calculations 


STEADY  FLOW  IN  THE  LAMINAR  LAYER  BETWEEN 
TWO  PARALLEL  STREAMS 


The  basic  steady  flow  -under  discussion  is  a  boundary-layer  flow  with 
no  body  forces  and  no  pressure  gradient.  The  basic  equations  are  (see 
list  of  symbols  in  the  appendix) : 


p*u*  +  p*v*  Su* 
8x* 


8v* 


-U*  ^ 

8y*\  8y*/ 


(l) 


— (p*u*)  +  iL(p*v*)  =0  (2) 

8x*  By* 


*  * 
P  u 


*  *  8  /  m*\ 

+  P  V  8y*  'Cp  ' 


=  8  L*  8t*\ 

8y*  \  8y*/ 


(3) 
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the  pressure  "being  a  constant  throughout  the  field.  In  the  case  of  a 
homogeneous  incompressible  fluid,  with  the  streams  at  the  same  tempera¬ 
ture,  the  temperature  may  he  taken  as  constant  throughout  the  field, 
n.nrl  the  first  two  equations  can  he  solved  for  the  velocity  distribu¬ 
tions  u*(x*,y*)  and  v*(x*,y*) .  In  the  case  of  a  compressible  gas, 
the  integration  has  to  be  carried  out  for  individual  cases.  However, 
if  the  Prandtl  number  Cpp*/k*  is  equal  to  unity,  it  is  known  that 

there  is  a  quadratic  relation  of  the  type 


CriT*  +  i  u*2  =  A  +  Bu* 


(4) 


between  the  temperature  and  the  velocity,  and  one  is  again  essentially 
dealing  with  two  distributions  u*(x*,y*)  and  v*(x*,y*).  Indeed,  the 
constants  A  and  B  are  given  in  terms  of  the  conditions  in  the. 
parallel  streams  as  follows: 


"1 


k 


(5) 


A  =  -  |  UjU2  +  c^TjU.,  -  T2Ul)/(u2  -  Ux) 


If  the  total  enthalpy  in  the  two  streams  1  and  2  is  the  same,  that  is, 
if 


cpTl  +  |  Ul2 


=  cpT2 


+  iUc 


(6) 


then  B  =  0,  and  the  total  enthalpy  is  constant  throughout  the  whole  field: 


cpT*  +  |  u*2  =  C 


(7) 
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This  is  a  particularly  simple  case,  to  which,  however,  all  other  cases 
can  he  easily  reduced.  This  is  done  by  rewriting  equation  (4)  into 
the  form 


cpT*  +  i(u*  -  B)2  =  A  +  B2/2 
2 


One  need  only  consider  an  observer  moving  with  the  speed  B  and  con- 
sider  the  relative  velocity  u*  -  B.  Thus,  the  solution  for  general 
specified  values  of  Up,Tp,  and  U2,T2  can  be  derived  from  the  iso- 

energetic  solution  with  boundary  conditions  Up  -  B,Tp  and  U2  -  B,T2 
by  simply  adding  the  constant  B  to  the  u*-component  of  the  velocity. 1 

In  the  following  discussions,  isoenergetic  basic  solutions  will  be 
referred  to  often;  however,  it  should  be  kept  in  mind  that  by  the  con¬ 
sideration  of  a  moving  observer  the  general  case  may  be  obtained.'  This 
transformation  is  not  restricted  to  the  steady  flow  but  applies  to  the 
consideration  of  the  disturbances  as  well.  Thus,  if  all  the  cases  of 
constant  total  enthalpy  are  calculated,  all  the  other  cases  are  also 
known . 

So  far,  the  viscosity  coefficient  may  depend  on  the  absolute  tem¬ 
perature  in  any  manner.  If  there  is  direct  proportionality  of  these 
quantities,  the  solution  in  the  compressible  case  can  be  expressed  in 
terms  of  that  in  the  incompressible  case.  These  relations  are  well- 
known  and,  in  the  following  discussion,  only  the  results  relevant  to 
this  case  will  be  given. 


Incompressible  Case 
For  the  incompressible  case: 


U*  =  Upf  1  (  T]  ) 


V*  =  Upjj-  ^  Vp  |UpX*  (qf  ’ 


lit  is  to  be  noticed  that  in  certain  cases  with  T2/Tp  <1  the 
boundary  value  U2  -  B  for  the  corresponding  isoenergetic  problem  may 
be.  negative  while  Up  -  B  is  still  positive.  Thus  the  isoenergetic 
problem  may  not  be  physically  significant.  However,  for  the  purposes 
of  theoretical  analysis,  this  point  is  not  important. 
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where 


and  f (tj )  satisfies  the  differential  equation 


(10) 


ff ' '  +  2f  ”  '  =  0 


(ID 


with  the  conditions 


f  '(tj)  — >  1  as  T) 


f -(tj)  — >  ^2  as  Tj  > 
U1 


(12) 


A  third  condition  is  arbitrary  up  to  a  translation  along  the  q-axis . 

For  a  typical  scale,  the  momentum-boundary-layer  thickness  9*  may 
be  introduced,  which  is  defined  by 


foo 

p*(ul 


(13) 


Then 


(14) 


for  an  incompressible  fluid.  It  will  be  seen  later  that  the  same  for¬ 
mulas  apply  for  a  compressible  fluid  in  the  isoenergetic  case. 

Calculations  made  by  Gortler  (reference  8)  for  the  turbulent  mixing 
region  can  be  easily  adopted  for  the  purpose  at  hand.  There  is  only  a 
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slight  difference  in  the  method  of  representation.  To  convert  his 
function  F'(^)  to  the  notation  of  this  report,  the  following  rela¬ 
tions  should  he  used 


f '(n)  =  J£—  F'(g) 
1  +  X 

*1  =  2  Nl1  +  M 

where 


Ui  +u2 


(l5a) 


(15b) 


The  converted  results  are  given  in  table  I  and  figure  1. 

Compressible  Case 

For  the  compressible  case  with  constant  total  enthalpy,  p*  oc  T*, 

where 

u*  =  U1f>(0  (16) 

and 

T*  =  T±(l  +  2LlJ:  M^Jl  -  f'2(sj  (l7) 

with 

^1  -  21lJ:  Mx2[l  -  f,2(^}  d£  (18) 

where 


n  =  y*/Jvix*/ui 
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The  momentum  thickness  0*  can  he  defined  by  the  same  equation  (equa¬ 
tion  (13))  and  can  be  shown  to  have  the  same  value  (equation  ( 14) )  in 
dimensionless  form.  Velocity  and  temperature  distributions  in  the  case 
Ug  =  0  have  been  calculated  for  several  Mach  numbers;  these  are  tabu¬ 
lated  in  table  II  and  plotted  in  figure  2. 


A  GENERAL  RESTRICTION  FOR  EXISTENCE  OF 
NEUTRAL  SUBSONIC  DISTURBANCES 


Before  going  into  the  general  theory  of  stability,  a  preliminary 
discussion  of  the  stability  of  the  mixing  zone  will  be  given.  It  will 
be  found  below  that  the  results  thus  obtained  give  quite  an  adequate 
description  of  the  general  dependence  of  the  stability  of  the  mixing 
zone  on  the  Mach  numbers  of  the  streams.  It  will  be  shown  that,  if 
the  average  Mach  number  "M  of  the  relative  motion  (defined  by  equa¬ 
tions  (25)  and  (26))  exceeds  the  value  2,  subsonic  disturbances  in  the 
sense  of  reference  1  cannot  exist,  and  the  mixing  zone  may  be  expected 
to  be  stable. 

Consider  two  parallel  streams  at  speeds  Up  and  U2,  tempera¬ 
tures  Tp  and  Tg,  and  Mach  numbers  Mp  and  Mg-  For  definiteness, 
take  Up  -  Ug  >  0.  Consider  an  observer  moving  with  the  speed  Ug. 
Then  the  streams  appear  to  have  speeds  Up  -  Ug  and  0,  while  their 
temperatures  are  obviously  not  changed.  If 


c'  =  c  -  Ug 


(19) 


denotes  the  speed  (relative  to  the  moving  observer)  of  a  wavy  motion 
propagating  in  the  direction  of  the  stream  (c'  positive  or  negative), 
then  the  conditions  for  subsonic  disturbances  are 


(Up  -  Ug)  -  c  <  ap 


(20) 


| c  '  |  <  a2  (21) 

where  ap  and  a2  are  the  velocities  of  sound  in  the  streams.  It  is 
clear  that  subsonic  disturbances  can  always  exist  if  Up  -  Ug  is  less 
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than  a-j_.  For  -  Ug  >  a^,  condition  (20)  can  he  satisfied  only 
with  c '  >  0 .  Thus 

0  <  c'  <  a2,  that  is,  0  <  e  -  U2  <  a2  (22) 

Now,  note  that  condition  (20)  is  actually  , 

|U1  “  c|  <  al 

and  hence 


-a1  <  %  -  c  <  ax  (23) 

Adding  relations  (22)  and  (23),  it  is  found  that 


U-J_  —  Ug  < 


(24) 


is  a  necessary  condition  for  the  existence  of  subsonic  disturbances. 
By  introducing  the  average  velocity  of  sound 


a  —  —  ^a^  +  a2^  (25) 

and  the  average  Mach  number  of  relative  motion 

'  -  < 

M  =  (u-l  -  U2)|a  (26) 


defined  with  respect  to  this  average  sound  velocity,  it  may  be  concluded 
that  subsonic  disturbances  cannot  exist  (and  the  mixing  zone  may  be 
expected  to  be  stable)  if  the  average  Mach  number  M  of  the  relative 
motion  exceeds  the  value  2,  as  stated  at  the  beginning  of  this  section. 
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In  terms  of  the  Mach  numbers  of  the  two  streams,  condition  (24) 
states  that  if 


«l  > 1  +  *  «a) 


(27) 


Subsonic  disturbances  cannot  exist,  and  the  motion  is  stable.  (See 

fig-  3-) 

In  the  case  of  constant  total  enthalpy,  the  temperature  ratio 
Tg/Ti  can  be  expressed  as 


t2 


1  +  (7  -  l)Mj2/2 

Tl  1  +  (7  -  l)M22/2 

Then  condition  (27)  becomes 

1  +  (7  -  1)M12/211/2 


Mx  >  1  + 


The  critical  condition  is 


1  +  (7  -  l)M22/2 


4  +  MZ) 


(28) 


(29) 


M1  =  1  + 


1  +  (7  -  l)M12y2 


1  +  (7  -  l)M22/2 


1/2 


(l  +  M2) 


Removing  the  square  root, 

1  +  (7  -  1)M12^2  -  l)2 


1  + 


(7  -  l)M22/2  (M2  +  L)‘ 


(30) 


This  form  suggests  an  obvious  solution,  Mx  =  -Mg. 
can  then  be  easily  obtained  as 


Another  solution 


Ml  = 


(3  -  7 )Mg  +  4 
(3  -  7)  -  2(7  -  1)M2 


(31) 
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By  substituting  in  equation  (30),  Mq  =  -Mg  is  found  to  be  an  extra¬ 
neous  solution.  The  only  solution  of  equation  (30)  is  then  given  by 
equation  (31). 

Thus  the  flow  is  stable  if 


(3  -  7) M2  +  4 

Mx  >  - 

(3  -  7)  -  2(7  -  l)Mg 


(32) 


The  curve  for  equation  (31)  is  shown  in  figure  4,  where  the  region  to 
the  left  of  curve  A  is  the  region  of  possible  instability.  Note  the 
symmetry  with  respect  to  the  line  Mq  =  -Mg,  corresponding  to  a  change 

of  the  positive  direction  of  the  x-axis  in  the  physical  problem.  The 
line  M^_  =  Mg  is  drawn  to  take  care  of  the  condition  Uq  -  Ug  >  0, 

since  Uq  =  Ug  when  M-^  =  Mg.  The  asymptotes  of  the  curve  for  equa¬ 
tion  (31)  are 

Mg  =  (3  -  7)/2(7  -  1)  (33) 


and 

Mq  =  -(3  -  7)/2(7  "  1)  (34) 

With  7  =  1.4,  these  values  are  ±2.  The  intercepts  are 


Mq  =  4/(3  -  7) 

and 

Mg  =  -4/(3  -  7) 


(35) 

(36) 


With  7  =  1.4,  these  values  are  +2.5.  (Cf.  reference  9-) 
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GENERAL  STUDY  OF  SMALL  DISTURBANCES  IN  A  NEARLY  PARALLEL 
FLOW  FIELD  IN  A  COMPRESSIBLE  FLUID 


The  general  theory  of  a  small  disturbance  in  a  field  of  nearly  par¬ 
allel  flow  of  a  gas  has  been  developed  in  reference  1.  The  rigorous 
mathematical  proof  and  improvement  of  the  theory  are  given  in  detail  in 
references  5  and  6.  In  this  report,  merely  the  main  conclusions  and 
their  physical  interpretations  are  outlined  without  going  into  the 
details .  Applications  of  the  theory  to  the  specific  case  at  hand  will 
be  discussed  in  some  detail. 

Consider  a  nearly  parallel  stream  with  dimensionless  velocity  and 
temperature  distributions  w(y)  and  T(y).  The  neglect  of  the  depend¬ 
ence  of  these  quantities  on  x  and  the  omission  of  the  y-component  of 
the  basic  flow  can  be  justified  by  detailed  investigations.  The  line¬ 
arized  differential  equations  for  small  disturbances  then  possess  solu¬ 
tions  of  the  type 


u'  =  Re[f(y)eia(x"ct)] 
v'  =  Re  ^(y)eia(x-cty| 

p'  =  Re  [£(y)e^a^x-c^^] 

p'  =  Ref(y)eia^x-Ct3 


T'  =  Re 


f(y)eia^x-ct3 


(37) 


where  u',  v',  p',  p',  and  T'  are  the  perturbations  of  the  two  com¬ 
ponents  of  velocity,  the  density,  the  pressure,  and  the  temperature,  all 
in  a  suitably  defined  dimensionless  form.  The  constants  a  and  c  are, 
respectively,  the  real  wave  number  and  the  complex  wave  speed. 

The  differential  equations  for  the  amplitude  functions  f (y) ,  0(y), 
r(y) ,  n(y) ,  and  @(y)  are  rather  complicated.  However,  if  the  effects 
of  viscosity  and  heat  conduction  are  neglected,  they  become  a  relatively 
simple  system,  which  can  be  reduced  to  the  following  single  differential 
equation  for  0(y): 


The  boundary  conditions  are  that  the  disturbance  should  be  bounded 
as  y  — >  ±oo. 


The  inviscid  system  would  have  given  a  well-defined  characteristic- 
value  problem  if  it  were  not  for  the  fact  that  differential  equation  (38) 
has  a  singularity  at  the  point  y  =  yc  where  w(y)  =  c.  This  singu¬ 
larity  disappears  only  when  d(pw')/dy  =  0  at  the  same  point.  Other¬ 
wise,  a  solution  of  equation'  (38)  has  a  logarithmic  singularity  at 
y  =  yc  and  the  characteristic-value  problem  associated  with  this  equa¬ 
tion  becomes  indeterminate  until  the  proper  branch  of  the  solution  is 
determined . 

The  determination  of  the  proper  branch  of  the  solution  and  its 
associated  physical  interpretation  is  one  of  the  most  delicate  points 
in  the  theory  of  hydrodynamic  stability.  The  mathematical  analysis  of 
the  solutions  of  the  complete  viscous  equations  and  their  limiting  solu¬ 
tions  will  be  made  first  before  discussing  their  physical  interpretation. 

The  complete  system  of  viscous  equations  can  be  shown  to  be  equiva¬ 
lent  to  a  system  of  six  linear  equations  of  the  first  order  in  six 
unknowns.  Thus,  there  are  six  independent  solutions.  These  solutions 
have  been  formally  obtained  as  asymptotic  series  in  reference  1,  and 
their  rigorous  mathematical  investigation  has  been  carried  out  in 
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references  5  and  6.  It  is  found  that  two  of  the  six  solutions  can  be 
expressed  in  asymptotic  series  of  the  form 


+  —  f(X)  +  . 
X2 


V 

0  =  0(0)  +  i_  0(1)  +  .  .  . 

X2 


(40) 


where  f(®),  0(<")),  .  .  .  are  the  inviscid  solutions  satisfying  equa¬ 

tions  (38)  and  (39)  and  X2  =  aR,  where  R  is  the  Reynolds  number  based 
on  the  thickness  of  the  mixing  region.  Thus,  the  formal  limit  of  equa¬ 
tions  (40)  does  approach  the  inviscid  solution,  but  a  complete  study  of 
these  equations  also  carries  the  knowledge  of  the  proper  branch  to  be 
used. 


Four  other  solutions  of  the  complete  system  of  viscous  equations 
are  of  the  form 


f  =  F  exp^XQ-jJ 
0=0  exp(XQi) 


i  =  1,  2,  3,  4  y 


(41) 


where 


>  (42) 


and  F  and  0  can  be  expressed  as  power  series  of  X,  involving  only 
a  finite  number  of  positive  powers. 
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In  the  establishment  of  these  asymptotic  solutions,  it  is  shown  that 
the  lines  ^e^Q-^j  =  0,  Re^Q^)  =  0,  and  Re^Q-^  -  =  0  are  of  interest. 

The  geometry  of  these  lines  relative  to  the  point  yc  (where  w  =  c)  and 

the  real  axis  of  the  y-plane  are  shown  in  figure  5.  There  are  asymptotic 
solutions  expressed  in  equations  (40)  and  (4l)  which  maintain  the  same 
analytical  expression  on  the  two  sides  of  the  dotted  lines.  However,  in 
crossing  the  solid  lines,  they  generally  change  their  behavior.  Thus, 
the  following  conclusion  may  be  drawn:  The  proper  branch  of  the  multiple¬ 
valued  asymptotic  solutions  is  obtained  by  taking  a  path  in  the  complex 
y-plane  below  the  point  y  =  yc  (in  this  case  w(y)  is  monotonically 
increasing  along  the  real  axis). 

This  is  the  branch  taken  in  reference  1.  By  examining  the  behavior 
of  solutions  of  the  type  given  in  equations  (4l),  it  can  be  easily  shown 
that  they  diverge  for  either  positive  or  negative  large  values  of  y. 

These  solutions  should,  therefore,  be  rejected  in  the  present  problem. 

The  effect  of  viscosity  is  then  to  be  obtained  through  solutions  of  the 
form  of  equations  (40) . 


GENERAL  STABILITY  CHARACTERISTICS  IN 
THE  INVISCID  CASE 


Much  of  the  discussion  of  the  inviscid  case  in  reference  1  applies 
to  the  present  case.  However,  as  noted  above,  in  the  present  problem, 
the  possibility  can  be  more  readily  realized  that  subsonic  disturbances 
may  not  exist  at  all,  and  that  the  motion  may  then  be  expected  to  be 
completely  stable  with  respect  to  small  disturbances .  Another  main 
difference  lies  in  the  "steady"  deviation,  that  is,  solutions  of  equa¬ 
tion  (38)  with  a  =  0.  In  the  present  discussion,  the  general  line  of 
discussion  in  reference  1  will  be  followed.  However,  the  difference  in 
the  boundary  conditions  often  causes  a  difference  in  the  method  of 
analysis.  The  arguments  are,  therefore,  presented  in  some  detail.  Also, 
it  will  be  shown  that  disturbances  having  wave-lengths  slightly  longer 
than  that  of  the  neutral  subsonic  disturbances  of  finite  wave  length 
are  unstable  while  those  with  slightly  shorter  wave  lengths  are  the 
stable  ones.  This  conclusion  also  applies  to  the  case  of  the  boundary 
layer,  but  it  was  not  obtained  in  reference  1.  The  analysis  also  leads 
to  an  approximate  estimation  of  the  dependence  of  amplification  on  wave 
length.  This  will  be  used  for  the  calculation  of  the  amplification  of 
the  disturbances  in  the  section  "Self -Excited  Oscillations." 

Some  general  analytical  properties  of  the  solutions  of  equation  (38) 
will  be  first  summarized,  particularly  for  the  case  of  subsonic  disturb¬ 
ances  .  In  this  case, 
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>  0 


(43) 


for  "both  free  streams.  For  large  positive  values  of  y,  equation  (38) 
may  t>e  approximated  by 


0"  -  Pl20  =  0  with  P-j2 


(44) 


For  large  negative  values  of  y,  it  may  be  approximated  by 


2  2 

,  7u2  V 

L 

=  0  with  p2  =  a 

't2 

(45) 


Since  both  px2  and  Pg2  are  positive,  the  solution  0(y)  is  exponen¬ 
tial  in  nature  for  large  values  of  y  in  the  case  of  subsonic  disturb¬ 
ances,  with  the  exception  of  the  case  a  =  0.  In  that  case,  two  inde¬ 
pendent  solutions  of  equation  (38)  are 


=  w  -  c 


(46) 


and 


0 


2 


(w  - 


(47) 


The  first  solution  is  bounded  while  the  second  varies'  linearly  with  y 
for  large  values  of  y. 

For  any  value  of  c,  the  bounded  solution  (equation  (46))  corre¬ 
sponds  to  no  disturbance  at  infinity.  In  fact,  for  a  =  0,  the  disturb¬ 
ance  v’  is  identically  zero,  by  equation  (37).  The  other  components 
of  the  disturbance  are  given  by  equations  (39),  and  it  can  be  easily 
verified  that  they  all  vanish  at  large  distances. 

It  can  be  shown  from  the  general  nature  of  the  temperature  and 
velocity  distributions  that  condition  (43)  holds  throughout  the  mixing 
zone  if  it  holds  in  both  free  streams.  It  is  then  obvious  that 
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equation  (38)  has  a  singularity  only  at  w  -  c  =  0.  Actually,  even 
if  T  -  M-|_2(-w  -  e)^  =  0,  it  only  gives  rise  to  an  apparent  singularity. 

Analysis  of  the  solution  in  the  neighborhood  of  the  point  y  =  yc 
(where  w  =  c)  gives  the  following  two  solutions: 


0i  =  (y  -  yc)si(y  -  yc) 

02  =  s2(y  -  yc)  +  K0i  l0§e  (y  -  yc) 


(48) 


where 


T„ 


2  1—  /  n 


K  = 


(V)3 


d  /w1^ 


(49) 


and  g-j^  and  g2  are  power  series  in  (y  -  yc)  with  g1(0)  =  wc'  £  0 

and  g2(0)  =  Tc/wc'  ^  0.  The  proper  branch  of  the  logarithmic  function 

is  to  be  taken  in  accordance  with  the  method  discussed  in  the  preceding 
section. 

For  real  values  of  c,  it  can  be  shown  that  the  Reynolds  shear 
stress 


t  =  -pu'v' 


(50) 


is  a  constant  except  for  a  possible  jump  at  y  =  y  . 
jump  is 


In  fact,  this 


[x]  =|*k(vo’)2|0c|2/tc 


(51) 


when  one  passes  from  yc  -  0  to  yc  +  0.  Thus,  for  a  neutral  disturb¬ 
ance  with  a  ^  0,  the  condition 


K  =  0 


(52) 


must  be  satisfied,  since  t  =  0  for  y  — ±00  and,  therefore,  does  not 
have  any  jump. 
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Thus, 


for  the  existence  of  a 


subsonic  disturbance,  the  quantity. 


_d_/l  dw\ 
dy  \T  dy/ 

must  vanish  at  some  point  in  the  field;  furthermore,  the  corresponding 
value  wg  of  w  must  be  subsonic  relative  to  both  streams;  that  is. 


1  - 


I 


(54) 


The  above  reasoning  can  be  applied  to  the  case  of  the  boundary  layer; 
a  somewhat  different  argument  was  used  in  reference  1.  A  little'  calcu¬ 
lation  will  show  that,  in  general,  condition  (53)  will  be  satisfied, 
although  condition  (54)  may  not.  If  the  latter  is  also  satisfied,  then 
there  actually  exists  a  subsonic  disturbance  with  c  =  wg .  To  prove 
this,  equation  (38)  is  rewritten  in  the  form 


_d_ 

ay 


(55) 


where 


h(y)  =  [t  -  M1z(w  - 


q(y)  =  -±-  ±(h  &) 

w  -  C  dy\  dy/ 


(56) 


With  c  =  ws,  q(y)  is  regular  all  along  the  real  axis.  This  is  a 

2 

characteristic-value  problem  for  the  parameter  k  =  -a  and  is  associ¬ 
ated  with  the  variational  principle 


5 


dy  =  0 


(57) 
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with 


(58) 


O 

The  least  value  for  -a  is  then  given  "by  the  minimum  value  of  the 
ratio 


dy  =  Constant 


(59) 


for  all  functions  f(y)  such  that  the  integrals  in  this  ratio  are  con¬ 
vergent.  Now,  it  is  necessary  only  to  show  that  for  certain  functions 
the  ratio  of  the  integrals  is  negative. 

In  the  choice  of  such  a  function,  it  is  necessary  for  only  the 
numerator  in  equation  (59)  to  be  convergent.  The  convergence  of  the 
integral  in  the  denominator  is  immaterial,  for  one  can  always  modify 
the  function  at  sufficiently  large  values  of  |y|  so  that  the  denomi¬ 
nator  becomes  convergent  without  altering  the  sign  of  the  numerator. 

To  choose  a  function  so  that  the  integral 


1=1  (hf '2  +  qf2)  dy 


(60) 


is  negative,  it  is  first  noted  that 


1=0  for  f  =  w  -  c 


(6l) 


This  follows  from  the  fact  that  f)  =  w  -  c  is  a  solution  of  equation  (55) 
when  a  =  0.  It  can  also  be  directly  verified.  Obviously,  the  value  of 
I  is  not  changed  if  f  is  now  changed  to  | w  -  c | .  A  further  modifi¬ 
cation  of  the  function  would  yield  the  desired  result.  In  the  neighbor¬ 
hood  of  yc,  |  w  -  c  |  is  small,  but  w'  is  finite.  If  in  this  neighbor¬ 
hood  the  function  f(y)  =  |w  -  cj  is  replaced  by  a  horizontal  straight 
line,  the  integrand  in  equation  (60)  is  certainly  decreased  (since  |  >0 
and  a  change  in  f '  causes  a  larger  change  in  the  integral  than  that 
in  f),  and  the  resultant  integral  is  negative.  This  completes  the  proof 
desired. 
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Thus,  the  sufficiency  of  conditions  (53)  and  (54)  for  the  existence 
of  subsonic  disturbances  is  also  established.  Let  the  corresponding 
value  of  a  be  denoted  by  as . 

Next,  it  will  be  shown  that  for  a  slightly  less  than  as  the 
disturbances  are  unstable.  To  prove  this  the  quantity  (dc/dk)g  will 
be  calculated,  and  its  imaginary  part  will  be  shown  to  be  positive. 

Consider  a  characteristic  function  0(y;k,c)  of  equation  (55) 
corresponding  to  a  given  value  of  k.  As  c  changes,  k  also  changes. 
If  equation  (55)  is  differentiated  with  respect  to  k,  the  following 
equation  is  obtained 


dy  \  dy 


where 

juM  +  MS£  (63) 

8k  8c  dk 

Now,  multiply  equation  (62)  by  0,  equation  (55)  by  0,  subtract  the 
results,  and  then  integrate  with  respect  to  y  along  a  path  in  the  com¬ 
plex  plane  which  leads  from  y  =  -00  on  the  real  axis  to  y  =  00  on  the 
real  axis  but  passes  below  the  point  yc.  The  following  equation  is  then 
obtained  for  dc/dk: 


dc  T  8h/d0\2  +  8q  jZ 
dk  J_ro  [8c\dy/  8c  _ 

So  far  equation  (64)  is  general.  Now,  specialize  to  the  case  of 
the  neutral  disturbance  in  question.  Most  of  the  integration  can  then 
be  carried  out  along  the  real  axis,  with  real  resultant  values.  However, 
the  integral 


dy  = 


T 


dy 


(64) 


dy 


(65) 
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must  be  evaluated  along  an  indented  path  with  a  small  detour  below  yc , 

since  6q/8c  has  a  pole  at  that  point.  Calculations  show  that  the 
imaginary  part  of  J  is 


if  this  quantity  does  not  vanish.  Indeed,  this  quantity  can  be  seen 
to  be  negative.  Thus,  equation  (64)  yields  a  relation  of  the  form 


—  (P  +  iQ)  =  R  >  0 
dk 

where  Q  =  Im(j) .  Thus, 


dc 

dk 


R 

P2  +  q2 


(P  -  iQ) 


(67) 


has  a  positive  imaginary  part,  as  required. 

In  the  incompressible  case,  equation  (64)  reduces  to 


dc 

dk 


(68) 


and  the  imaginary  part  of  dk/dc  is 


Im(dk/dc) 


(69) 
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NEUTRAL  OSCILLATIONS 


From  the  discussions  of  the  last  section,  another  limitation  is 
imposed  on  the  occurrence  of  neutral  subsonic  disturbances .  Since  the 
wave  speed  of  the  neutral  disturbance  must  be  equal  to  the  flow  speed 
ws  at  the  point  where 

dy  V  dy/ 

subsonic  disturbances  cannot  occur  if  the  speed  ws  is  supersonic  rela¬ 
tive  to  either  stream.  Calculations  of  ws  are  made  for  the  case  of 
uniform  total  enthalpy.  Instead  of  using  equation  (70)  directly,  it 
is  found  convenient  to  transform  it  into  the  form 


-1 

=  0  with  f  =  f(£)  and  w  =  f'(^)  (71) 


where  i|r  is  the  stream  function,  since  the  relation  w(i|r)  .is  identical 
with  that  in  the  incompressible  case.  The  results  of  these  calculations 
are  shown  in  table  III.  In  comparing  these  values  against  condition  (54) 
it  is  found  that  the  condition 


8^ 

a? 


1  + 


(7  - 


is  never  violated,  but  the  values  of  ws  below  the  solid  horizontal 
lines  in  the  table  are  too  high  to  satisfy  the  condition 


Thus,  for  a  given  Mach  number  in  the  stream  Mg,  the  Mach  number  M-j_ 
in  the  other  stream  can  become  so  high  that  the  speed  of  flow  corre¬ 
sponding  to  condition  (70)  becomes  supersonic  relative  to  the  slower 
free  stream.  Subsonic  disturbances  then  cannot  exist. 
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The  critical  case  of  sonic  disturbance  is  reached  when 


Calculations  for  this  case  are  tabulated  in  table  IV.  The  corresponding 
values  of  Mp  and  M2  are  plotted  as  curve  B  in  figure  4  to  mark  the 
limit  of  stability.  Only  the  region  to  the  left  of  this  curve  can  have 
subsonic  disturbances.  It  is  seen  that  the  condition  is  more  restrictive 
than  that  obtained  from  general  considerations  alone.  In  particular,  in 
the  case  where  one  stream  is  at  rest,  the  curve  B  shows  that  the  flow  . 
becomes  completely  stable  at  a  free-stream  Mach  number  of  1.7,  in  con¬ 
trast  to  the  value  2.5  given  by  curve  A.  Since  they  are,  however,  not 
very  much  different,  it  may  be  surmised  that  the  exact  distributions  of 
temperature  and  velocity  in  the  mixing  region  may  not  be  too  important 
in  determining  the  over-all  stability  characteristics. 

It  may  be  recalled  that  the  restriction  to  the  case  of  constant 
total  enthalpy  can  be  immediately  removed  by  considering  a  moving  obser¬ 
ver.  Thus,,  for  all  cases  with  Prandtl' number  equal  to  unity,  it  is  nec¬ 
essary  only  to  convert  the  values  of  the  speeds  Up,  JJg,  and  ws 

involved.  Condition  (70)  is  not  modified  by  the  reference  to  a  moving 
observer. , 

With  the  neutral  wave  speed  thus  determined,  equation  (55)  can  be 
integrated  to  give  the  amplitude  of  the  oscillations.  For  this  purpose, 
it  is  necessary  to  find  the  proper  value  for  a.  This  can  be  done  by 
several  trials,  with  the  first  approximation  given  by  the  ratio  of  the 
integrals  in  equation  (59) .  The  ratio  will  yield  the  characteristic 
value  only  when  the  function  f(y)  is  the  characteristic  function,  but 
it  is  known  that  any  reasonable  approximation  to  it  will  give  a  very 
close  approximation  to  the  characteristic  value. 

Calculations  of  these  neutral  oscillations  are  carried  out  for 
several  Mach  numbers  in  one  stream  with  the  other  at  rest.  These  are 
given  in  table  V  and  figure  6. 


SELF-EXCITED  OSCILLATIONS 


The  formulas  given  in  the  section  "General  Stability  Characteristics 
in  the  Inviscid  Case"  have  been  used  to  calculate  the  characteristics  of 
self-excited  oscillations.  It  is  found  that 


dc 

da 


=  0.177  -  0.209i 
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for  the  case  Mp  =  1  and 

—  =  0.093  -  0.287i 
do, 

for  the  incompressible  case.  The  results  are  shown  diagrammatically  in 
figure  7.  It  is  seen  that  the  extent  of  amplification  is  fairly  large. 
There  is  also  indicated  a  decrease  of  amplification  with  increasing  Mach 
number . 

It  would  be  easy  to  calculate  the  amplitude  function  of  these  self- 
excited  oscillations  by  using  the  characteristic  values  obtained  above. 
This  was  not  carried  out  because  of  limitations  in  time. 


NATURE  OF  OSCILLATIONS  IN  LIMIT  OF 
INFINITE  REYNOLDS  NUMBER 


Calculations  in  the  section  "Neutral  Oscillations"  indicate  that 
the  neutral  oscillation  has  quite  a  simple  amplitude  distribution.  In 
fact,  it  does  not  show  any  node.  The  amplified  disturbances  are  expected 
to  show  similar  characteristics.  The  mathematical  theory  (references  5 
and  6),  however,  indicates  that  the  damped  oscillations  behave  in  a  much 
more  complicated  manner.  It  is  concluded  that  there  is  always  a  finite 
viscous  region  in  the  interior  of  the  fluid,  no  matter  how  large  the 
Reynolds  number  may  be.  In  fact,  a  minimum  width  of  this  region  is 
determined.  This  has  to  do  with  the  crossing  of  the  solid  lines  by  the 
real  axis  of  y  in  figure  5.  The  solution  in  the  part  of  the  real  axis 
between  the  solid  lines  shows  exponential  behavior  -  and  therefore  viscous 
nature  -  if  the  solution  in  the  outside  parts  shows  the  inviscid  behavior. 
This  type  of  conclusion  has  been  reached  in  reference  1.  However,  it 
was  possible  only  to  suggest  that  such  viscous  behavior  would  occur  at 
the  solid  lines.  The  improved  theory  shows  that  it  must  occur  through¬ 
out  the  region  in  between.  It  is  important  to  note  that  the  coirplex  con¬ 
jugate  of  an  amplified  solution  does  not  represent  a  damped  oscillation 
and  vice  versa,  although  this  conclusion  can  be  easily  reached  by  a 
cursory  examination  of  the  inviscid  equations.  The  damped  oscillations 
do  not  satisfy  the  inviscid  equation  all  along  the  real  axis;  otherwise, 
.they  could  not  take  on  the  proper  branch  of  the  logarithm  as  specified 
in  the  section  "General  Study  of  Small  Disturbances  in  a  Nearly  Parallel 
Flow  Field  in  a  Compressible  Fluid."  They  exhibit  a  behavior  very  much 
like  that  of  the  vorticity  distribution  in  fully  developed  turbulence 
flow.  For  large  Reynolds  numbers,  there  is  one  part  of  the  space  where 
the  vorticity  is  highly  concentrated;  in  another  part,  there  is  very 
little  vorticity.  This  illustrates  the  two  kinds  of  limiting  behavior 
of  a  viscous  fluid  in  the  limit  of  infinite  Reynolds  number:  In  one 
part  of  the  field,  the  inviscid  behavior  is  approached;  in  another  part, 
it  becomes  highly  oscillatory  spatially. 
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CONCLUSIONS 


From  a  study  of  the  stability  of  the  mixing  of  two  parallel  streams 
in  a  gas,  the  following  conclusions  may  be  drawn: 

1.  If  the  relative  speed  of  the  two  parallel  streams  exceeds  the  sum 
of  their  velocities  of  sound,  subsonic  oscillations  cannot  occur,  and 
the  mixing  region  may  be  expected  to  be  stable  with  respect  to  small 
disturbances . 

2.  A  further  necessary  condition  for  the  possible  occurrence  of 
small  subsonic  disturbances  is  that  somewhere  in  the  field 


where  y  is  a  positional  coordinate  across  the  stream,  p  is  the  den¬ 
sity  of  the  gas,  and  w  is  a  dimensionless  velocity  distribution.  This 
condition  is  usually  satisfied  for  the  present  class  of  problems. 

3.  If  the  speed  of  the  flow  at  such  a  point  is  denoted  by  ws, 
then  the  field  of  flow  can  execute  a  neutral  wavy  oscillation  having 
a  finite  wave  length  and  propagating  with  the  speed  .  c  =  ws  if  and 
only  if  ws  is  subsonic  relative  to  both  streams .  There  is  no  other 

possible  neutral  oscillation.  This  leads  to  a  more  strict  condition 
of  stability  than  that  given  by  conclusion  1. 

L.  Under  the  above  conditions,  the  field  of  flow  can  execute 
amplified  wavy  oscillations  having  wave  lengths  longer  than  that  of 
the  neutral  oscillation.  Oscillations  having  shorter  wave  lengths  are 
damped,  (This  specific  form  of  the  conclusion  was  not  obtained  for 
the  boundary  layer  at  a  solid  surface  in  NACA  TN  1115, ,but  its  validity 
can  be  shown  by  the  present  method.)  The  extent  of  amplification  in 
such  cases  is  fairly  large. 

5.  At  large  Reynolds  numbers,  the  amplified  disturbances  are 
essentially  free  from  the  effect  of  viscosity.  On  the  other  hand, 
disturbances  with  finite  damping  are  expected  to  exhibit  a  highly 
oscillatory  behavior  over  a  finite  region  in  the  field  of  flow.  This 
is  similar  to  the  structure  of  the  vorticity  field  in  fully  developed 
turbulence . 

6.  For  the  case  of  constant  enthalpy  with  one  stream  at  rest,  the 
wave  length  of  the  neutral  disturbances  increases  with  increasing  Mach 
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number  in  the  other  stream.  The  flow  becomes  completely  stable  at  a 
free-stream  Mach  number  of  1.7*  This  is  more  restrictive  than  the 
value  2.5  obtained  by  applying  conclusion  1  to  the  present  case. 


Massachusetts  Institute  of  Technology 
Cambridge,  Mass.,  July  30>  1952 
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APPENDIX 

SYMBOLS 


The  quantities  hearing  subscript  1  in  the  last  column  are  the 
dimensional  quantities  in  the  first  free  stream.  Corresponding  quan¬ 
tities  in  the  second  stream  hear  a  subscript  2.  The  quantities  with¬ 
out  a  prime  satisfy  the  equations  of  steady  motion;  those  with  primes 
satisfy  the  disturbance  equations. 


Dimensional 

quantities 

Positional  coordinates: 
x* 


Dimensionless 

quantities 


Reference 

quantities 


x  l  =  2^2Vqx*/Uq 

y*  y  =  n/2\|2~  l 

Time: 

t*  t  1/  Uq 

Velocity  components  in  directions  of  x-  and  y-axes,  respectively: 

(y)+f(y)ela(x-ct)  Uq 

-a0(y)e^a'^x_c^)  Uq 


u*  +  U*’ 


w{ 


V*  +  V* 1 

Density  of  gas: 

p*  4.  p*  * 

Pressure  of  gas: 

p*  +  p*' 

Temperature  of  gas: 

ip*  t 


p(y)+*(y)eia(x-ct) 


T(y)+0(y)eia^x-ct^ 


Pi 


Pi 


Coefficient  of  viscosity  of  gas: 
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Dimensional  Dimensionless 

quantities  quantities 

Thermal  conductivity: 

k*  n/Pr 

Wave  number  of  disturbance: 

a*  =  2rt/x*  a  =  2jt/x 

Phase  velocity  of  disturbance: 

c*  c 

Specific  heat  at  constant  volume: 

cv  1 

Specific  heat  at  constant  pressure: 

CP  y 

Gas  constant  per  gram: 

R*  7-1 

Reynolds  number 


Reference 

quantities 


Cpdl 


Ul 


cv 


Cv 


R  =  PiUpl/pp 


Mach  number 


Prandtl  number 
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TABLE  I 

VELOCITY  DISTRIBUTIONS  IN  INCOMPRESSIBLE  CASE 

[For  definition  of  symbols,  see  equations  (10), 
(l4),  and  (l5b);  v  =  u*/u^j 


(a)  For  X  =  0.2  and  9  =  O.O98 


n 

V 

n 

V 

V 

-9.0 

O.667 

-2.6 

0.684 

1.4 

0.941 

-8.8 

-2.4 

.689 

1.6 

.952 

-8.4 

-2.2 

.695 

1.8 

.961 

-8.0 

-2.0 

.701 

2.0 

.969 

-7.6 

-1.8 

.709 

2.2 

.976 

-7.2 

-1.6 

.719 

2.4 

.982 

-6.8 

-1.4 

.729 

2.6 

.986 

-6.4 

-1.2 

.741 

2.8 

.990 

-6.0 

-1.0 

.754 

3.0 

•  993 

-5.6 

-.8 

.768 

3-2 

•  995 

-5.2 

.668 

-  .6 

.784 

3^ 

.996 

-4.8 

.668 

-.4 

.800 

3-6 

.997 

-4.4 

.668 

-.2 

.816 

3.8 

.998 

-4.0 

.669 

0 

.833 

4.0 

.999 

-3-8 

.670 

.2 

.851 

4.2 

.999 

-3-6 

.671 

.4 

.868 

4.4 

•  999 

-3-4 

.673 

.6 

.884 

4.6 

1.000 

-3-2 

.675 

.8 

.900 

4.8 

1.000 

-3.0 

.677 

1.0 

.915 

5.0 

1.000 

-2.8 

.680 

1.2 

.928 

32 
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TABLE  I 

VELOCITY  DISTRIBUTIONS  IN  INCOMPRESSIBLE  CASE  -  Continued 


(c)  For  X  =  0.6  and  9  =  0.573 


n 

w 

4 

w 

4 

H 

w 

-9.0 

0.250 

-2.6 

0.319 

1.4 

0.848 

-8.8 

.250 

-2.4 

•  331 

1.6 

.873 

-8.4 

.250 

-2.2 

.345 

1.8 

.896 

-8.0 

.250 

-2.0 

.361 

2.0 

.915 

-7.6 

.250 

-1.8 

•  379 

2.2 

•  933 

-7-2 

.250 

-1.6 

.398 

2.4 

•  947 

-6.8 

.251 

-1.4 

.420 

2.6 

•  959 

-6.4 

.251 

-1.2 

.444 

2.8 

.969 

-6.0 

.252 

-1.0 

.470 

3-0 

•  977 

-5.6 

.253 

-.8 

.498 

3-2 

.983 

-5.2 

.255 

-  .6 

.527 

3-4 

.988 

-4.8 

.258 

-.4 

.559 

3-6 

•  991 

-4.4 

.263 

-.2 

.591 

3.8 

•  994 

-4.0 

•  2.70 

0 

.625 

4.0 

•  996 

-3-8 

.274 

.2 

.659 

4.2 

•  997 

-3.6 

•  279 

.4 

.693 

4.4 

•  998 

-3-4 

.285 

.6 

•  727 

4.6 

•  999 

-3.2 

.291 

.8 

.760 

4.8 

•  999 

-3-0 

•  299 

1.0 

.791 

5.0 

1.000 

-2.8 

•  308 

1.2 

.821 

3^ 
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TABLE  I 

VELOCITY  DISTRIBUTIONS  IN  INCOMPRESSIBLE  CASE  -  Concluded 


(e)  For  X.  =  1.0  and  9  =  1.160 


n 

w 

n 

V 

n 

V 

-9.0 

0 

-2.6 

0.127 

1.4 

O.783 

-8.8 

.000 

-2.4 

.143 

1.6 

.816 

-8.4 

.001 

-2.2 

.161 

1.8 

.847 

-8.0 

.001 

-2.0 

.181 

2.0 

.874 

-7.6 

.002 

-1.8 

.202 

2.2 

.898 

-7.2 

.004 

-1.6 

.227 

2.4 

.919 

-6.8 

.006 

-1.4 

.253 

2.6 

•  937 

-6.4 

.009 

-1.2 

.282 

2.8 

.951 

-6.0 

.013 

-1.0 

•313 

3-0 

.963 

-5.6 

.018 

-.8 

.347 

3.2 

.973 

-5.2 

.024 

-.6 

.383 

3-4 

.980 

-4.8 

.033 

-.4 

.420 

3.6 

.986 

-4.4 

.042 

-.2 

.459 

3-8 

•  990 

-4.0 

.055 

0 

.500 

4.0 

•  994 

-3.8 

.062 

.2 

.542 

4.2 

•  996 

-3.6 

.070 

.4 

.584 

4.4 

.997 

-3.4 

.079 

.6 

.626 

4.6 

■  998 

-3.2 

.089 

.8 

•.668 

4.8 

•  999 

-3.0 

.100 

1.0 

.708 

5.0 

1.000 

-2.8 

.113 

1.2 

•747 
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TABLE  II 

VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS 
IN  COMPRESSIBLE  CASE 


^One  stream  at  rest;  A,  =  1;  0  =  I.I60J 


(a)  M]_  =  0.5 


4 

w 

T 

■■ 

w 

T 

-12.6 

-12.2 

-11.8 

-11.4 

-11.0 

-10.6 

-10.2 

-9.8 

0 

1.050 

-1.4 
-1.2 
-1.0 
-.8 
-  .6 
-.4 
-.2 

0 

0.261 

.289 

.320 

.352 

.386 

.423 

.461 

.500 

1.047 

1.046 

1.045  • 

1.044 

1.043 

1.041 

1-039 

1.038 

-9A 

0 

1.050 

.2 

.540 

1.035 

-9.0 

.001 

.4 

.581 

I.033 

-8.6 

.001 

.6 

.622 

I.031 

-8.2 

.002 

.8 

.662 

1.028 

-7.8 

.003 

1.0 

.702 

1.025 

-7A 

.004 

1.2 

.740 

1.023 

-7.0 

.006 

1.4 

.776 

1.020 

-6 .6 

.009 

.  1.6 

.809 

1.017 

-6.2 

•  013 

1.8 

.840 

1.015 

-5.8 

.019 

2.0 

.868 

1.012 

-5.4 

.025 

2.2 

.892 

1.010 

-5.0 

•  033 

1.049 

2.4 

.914 

1.008 

-4.6 

.043 

2.6 

.932 

1.007 

-4.2 

.054 

2.8 

.9^7 

1.005 

-3.8 

.069 

3.0 

.960 

1.004 

-3.4 

.086 

3.4 

.978 

1.002 

-3.0 

.108 

3-8 

.989 

1.001 

-2.8 

.121 

4.2 

•  995 

1.001 

-2.6 

.136 

4.6 

.998 

1.000 

-2.4 

.152 

5.0 

•  999 

1.000 

-2.2 

.170 

5.4 

1.000 

1.000 

-2.0 

.190 

1.048 

5.6 

1.000 

1.000 

-1.8 

-1.6 

.211 

.235 

1.048 

1.047 

6.0 

1.000 

1.000 
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TABLE  II 


VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS 


IN  COMPRESSIBLE  CASE  -  Concluded 


(c)  Mj_ 


1.450 


1.449 

1.449 

1.449 

1.448 

1.448 

1.447 

1.445 

1.444 

1.441 

1.438 

1.435 

1.433 

1.430 

1.426 

1.422 

1.418 

1.413 


O.3IO 

•334 

.359 

.385 

.412 

.440 

.470 

.500 

.531 

.564 

.596 

.629 

.662 

.694 

.727 

.758 

.788 

.817 

.844 

.869 

.891 

.911 

.929 

.957 

.976 

.988 

.994 

•  998 

•  999 
1.000 
1.000 


1.407 

1.400 

1.392 

1.384. 

1.374 

1.363 

1.351 

1.338 

1.323 

1.307 

1.290 

1.272 

1.253 

1.233 

1.213 

1.192 

1.171 

1.150 

1.130 

1.111 

1.093 

1.076 

1.062 

I.O38 

1.021 

1.011 

1.005 

1.002 

1.001 

1.000 

1.000 
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TABLE  III 

WAVE  SPEEDS  OF  POSSIBLE  NEUTRAL  SUBSONIC  DISTURBANCES 
^Values  of  ws  =  u*/U;j_  defined  by  equation  ( 70) J 


M1 

X 

0.2 

0.4 

0.6 

0.8 

1.0 

0 

O.83L 

0.732 

0.657 

0.613 

0.576 

.5 

.845 

•743 

.678 

.632 

.594 

1 

.851 

.762 

.704 

.664 

.633 

2 

.878 

.823 

.796 

.772 

.762 

5 

•  949 

.940 

.941 

.938 

.938 

10 

.984 

.981 

.981 

.982 

.982 

^5^ 


'  TABLE  TV 

CONDITIONS  ASSOCIATED  WITH  NEUTRAL  SONIC  DISTURBANCE 


Conditions 

X. 

• 

0.2 

0.4 

0.6 

0.8 

1.0 

Mi 

>10 

4.0 

2.7 

2.0 

1.7 

m2 

>1.92 

.90 

.44 

.167 

0 

c 

.91 

.82 

.78 

.73 
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TABLE  V 

WAVE  NUMBER  AND  AMPLITUDE  OF  NEUTRAL  OSC ILLATIONS 

[y  =  n/2^2,  e  =  1-160] 

(a)  Wave  number  of  neutral  oscillations 
for  various  Mach  numbers 


Mx 

Wave  number,  a 6 

0 

0.459 

.5 

.454 

1.0 

•  374 

1.5 

.324 
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TABLE  V 

WAVE  NUMBER  AND  AMPLITUDE  OF  NEUTRAL  OSCILLATIONS  -  Continued 


(b)  Amplitude  of  neutral  oscillations 
for  M^  =  0 


y 

0 

d0/dy 

y 

0 

d0/dy 

l  y 

0 

<10/ dy 

-6.4 

0.001 

0.001 

-1.0 

0.297 

0.269 

0.107 

-0.107 

-6.2 

.001 

.001 

-.9 

.325 

.281 

.087 

-.087 

-6.0 

.001 

.001 

-.8 

.354 

.291 

.071 

-.071 

-5.8 

.002 

.002 

-.7 

.383 

.296 

.058 

-.058 

-5 .6 

.002 

.002 

-.6 

.412 

.294 

3*0 

.048 

-.048 

-5.4 

.002 

.003 

-.5 

.442 

.286 

.039 

■  -,039 

-5.2, 

•  003 

•  003 

-.4 

.469 

.269 

.032 

-.032 

-5.0 

.004 

.004 

-•3 

.495 

.243 

.026 

-.026 

-4.8 

.005 

.005 

-.2 

.518 

.206 

.022 

-.022 

-4.6 

.006 

.007 

-.1 

.536 

.159 

.018 

-.018 

-4.4 

.008 

.008 

0 

.549 

.103 

4.2 

.014 

-.014 

-4.2 

.009 

.011 

.1 

.556 

.040 

4.4 

.012 

-.012 

-4.0 

.012 

.013 

.2 

.557 

-.028 

4.6 

.010 

-.010 

-3.8  • 

.015 

.016 

•3 

.550 

-.097 

4.8 

.008 

-.008 

-3-6 

.018 

.021 

.4 

.537 

-.162 

5.0 

.007 

-.007 

-3-4 

•  023 

.026 

.5 

.518 

-.220 

5.2 

.005 

-.005 

-3-2 

.029 

.032 

.6 

.494 

-.267 

5.4 

.004 

-.004 

-3-0 

.036 

.040 

.7 

.465 

-.303 

5.6 

.004 

-.004 

-2.8 

.045 

.050 

.8 

.434 

-.325 

5.8 

.003 

-.003 

-2.6 

.056 

.062 

•  9 

.401 

-•334 

6.0 

.002 

-.002 

-2.4 

.070 

•  077 

1.0 

.367 

-•332 

6.2 

.002 

-.002 

-2.2 

.087 

.095 

1.1 

.335 

-.322 

6.4 

.002 

-.002 

-2.0 

.108 

.116 

1.2 

■303 

-.305 

6.6 

.001 

-.001 

-1.9 

.120 

.129 

1-3 

.274 

-.283 

6.8 

.001 

-.001 

-1.8 

.134 

.142 

1.4 

.247 

-.260 

7.0 

.001 

-.001 

-1.7 

.149 

.156 

1.5 

.222 

-.237 

-1.6 

.165 

.172 

1.6 

•  199 

-.213 

-1.5 

.183 

.188 

1.7 

.179 

-.192 

-1.4 

.202 

.204 

1.8 

.161 

-.171 

-1.3 

.224 

.221 

1.9 

.145 

-.152 

-1.2 

.247 

.238 

2.0 

.131 

-.135 

-1.1 

.271 

.254 

2.1 

.118 

-.113 

(a)  Velocity  distributions,  w  =  U*/U]_. 


Figure  2.-  Velocity  and  temperature  distributions  in  compressible  case, 
y*  UX  -  Ug 

T]  =  - ;  A,  =  ^ — — rj- ■■  =  1.0;  My,  Mach  number  of  moving  stream. 

1  2 

M2  =  0. 
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* 


Figure  2.-  Concluded. 
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Figure  4.-  Stable  zones  for  isoenergetic  case.  7  =  1.40;  M]_  and  Mg, 

Mach  numbers  of  two  streams;  ws,  wave  speed  of  neutral  disturbance; 
Ui  and  U2,  velocities  of  two  streams. 
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(b)  For  =  0.5  and  cud  =  0.454. 

Figure  6.-  Amplitude  of  neutral  oscillations,  y  =  — i—  tj  =  — _ - _ • 

2  'fz  2V^2  /v1x*/u1 

amplitude;  cl,  wave  number;  0,  momentum-boundary-layer  thickness; 

0  =  I.160  in  t)-units . 


m 
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(b)  For  Mi  =  1. 


Figure  7--  Extent  of  amplification  at  infinite  Reynolds  number. 

7  1 


a,  wave  number;  wave  length  in  y-units;  y  =  1 


y* 


2/2  2/2 

R,,  Reynolds  number;  c,  dimensionless  complex  wave  speed  (cr  +  ic^). 
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